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Abstract 



The self-consistent chiral soliton of the Nambu- Jona-Lasinio model including the u, p 
and ai (axial-) vector meson fields besides the chiral angle is investigated. The resulting 
energy spectrum of the one particle Dirac Hamiltonian is strongly distorted leading to a 
polarized Dirac sea which carries the complete baryon number. This supports Witten's 
conjecture that baryons can be described as topological solitons. The exploration of the 
isoscalar mean squared radius of the nucleon exhibits that the repulsive character of the 
isoscalar vector field ui as well as the attractive features of the (axial-) vector mesons p 
and ai are reproduced in the Nambu- Jona-Lasinio model. The axial charge of the nucleon 
qa comes out far too small. This can be understood as an artifact of the proper time 
regularization prescription. 
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1. Introduction 

In the last decade a description of baryons as chiral solitons proved to be quite success- 
ful. The soliton picture of baryons is based on large Nc QCD considerations, N c being 
the number of colors. In the limit N c — > oo, QCD is equivalent to an effective theory 
of weakly interacting mesons Later Witten conjectured that in this effective meson 
theory baryons emerge as soliton solutions ||. Furthermore, in the low-energy limit this 
meson theory is dominated by the pseudoscalar mesons, the would-be Goldstone bosons 
of spontaneously broken chiral symmetry described in form of a non-linear cx-model. In 
order to implement the chiral anomaly this non-linear a- mo del has to be supplemented 
by the Wess-Zumino action. Introducing external gauge fields allows to extract the corre- 
sponding Noether currents. Especially, the baryon current arising from the Wess-Zumino 
term proves to be identical to the topological current thus supporting Skyrme's original 
work [Q. 

Incorporating vector and axial vector mesons in the non-linear cx-model the Skyrme 
model arises in the limit of infinitely heavy vector and axialvector mesons. In addition, 
the inclusion of vector and axial vector mesons with their physical masses cures several 
deficiencies of the Skyrme model, as e.g. incomplete description of electromagnetic prop- 
erties due to missing vector dominance {§] or wrong "high energy behavior" of tt — N phase 
shifts due to the higher order stabilization terms[[|]. This demonstrates the important role 
of vector and axialvector mesons in the soliton description of baryons. 

Although we have good confidence that the non-linear cx-model extended by vector 
mesons represents the low-energy approximation to the effective meson theory of QCD 
the actual effective theory is not explicitly known. Therefore Witten's conjecture cannot 
be proven within QCD without using further approximations. On the other hand, the 
low-energy form of the effective meson theory is almost entirely determined by chiral 
symmetry. Therefore the soliton picture of baryons should not depend on the details 
of chiral QCD dynamics. This suggests to investigate whether Witten's conjecture is 
fulfilled in simpler models for the quark flavor dynamics of QCD. For this purpose the 
Nambu-Jona-Lasinio (NJL) model is well suited. First of all, like QCD, it obeys chiral 
symmetry and, furthermore, it can be motivated as low-energy approximation to QCD||. 
Additionally, it can be bosonized by functional integral methods. The resulting effective 
meson theory is in satisfactory agreement with the low-energy meson data. In fact, its 
gradient expansion yields in leading order the linear cx-model and the Wess-Zumino action. 

The bosonized NJL model shows that the topological current arises in leading order 
gradient expansion from the vacuum part of the baryon current. The vacuum is defined 
with all negative fermion states occupied and the positive energy states being empty. 
Therefore the topological current can describe a non-zero baryon charge only if the vacuum 
is charged requiring that the valence quarks are bound into the Dirac sea. Only then the 
chiral field can carry a non-trivial baryon number. Thus the key assumption underlying 
the soliton picture of baryons is that the valence quarks have joined the Dirac sea. Within 
the NJL model we can test whether this assumption is fulfilled. 

In recent years the soliton of the NJL model has been extensively studied. First 
calculations were restricted to the chiral fieldJJ], || |J. Extending the model to also include 
the p meson provides only minor changes to the chiral soliton |K|. In both calculations the 
energy eigenvalue corresponding to the valence quarks is positive. The physical picture 
changes drastically if the chiral partner of the p meson, the ai axialvector meson, is added. 
Then the valence quarks become strongly bound and join the Dirac sea, i.e. the Skyrmion 



picture of baryons results JIT]. However, in these calculations the isoscalar vector meson 
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io was left out for technical reasons. The inclusion of the u meson should even favor the 
Skyrmion picture since the u introduces repulsion and hence increases the spatial size of 
the chiral field which in turn increases the binding of the quarks. 

The inclusion of the to meson leads to substantial technical complications because the 
uj field develops a non-zero time-like component in the chiral soliton. Due to the need for 
regularization the effective meson theory can be properly defined only via the continuation 
to Euclidean space. Thereby time-components of vector fields are continued to imaginary 
values, too. As a consequence the Euclidean action becomes complex and the evaluation 
of the NJL action for a soliton field leads to a non-Hermitean eigenvalue problem, see 
ref. for a proper treatment. In this calculation only the u meson and the chiral field 
have been included in a self-consistent soliton calculation. In the present paper we present 
the self-consistent soliton calculation in the NJL model with all low-lying two flavor vector 
and axial vector meson fields included. 



2. Bosonization of the NJL Model 

The starting point for the following considerations is the chirally invariant NJL model 



13, 14 



C = q(i$ - m°)q + 2g x J2 (j^y*?) 2 + (<?y*75<?) 



2 92 £ \ {q\l,qf + (gy757M?) 2 J , (2-1) 



wherein q denotes the quark spinors and m° the current quark mass matrix. Here we will 
work in the isospin limit, i.e. m° u = m° d = m°. The matrices A*/2 are the generators of the 



flavor group (A = y 2/iVyl). Furthermore we will restrict ourselves to two flavors (Nf = 
2) implying A* = t\ i = 0, . . . , 3. The coupling constants g± and g 2 will be determined 
from mesonic properties, cf. e.g. refs. |TJ], [T5|, [TB|, [T7| for the calculation of meson properties 
in the NJL model. 

Applying standard functional integral bosonization techniques the model Q2.1D can be 
rewritten in terms of composite meson fields M 

A F = Trlog(z^)=Trlog(z^ + f+ 7 54-(P fi S + P L st)), 

A m = [ d 4 x ( - J_tr(EtE - m°(S + st) + (m ) 2 ) + -Ltr^V + A^)) . 

(2.2) 

Here Pr^l = (1 ±7s)/2 are the projectors on right- and left-handed quark fields, respec- 
tively. Vfj, = J2l=oy^T a /2 and = J2l=o A®r a /2 denote the vector and axial vector 
meson fields. V® and A ^ are real in Minkowski space. The complex field S describes the 
scalar and pseudoscalar meson fields, SV,- = S a r°j/2 and Py = P a r° /2: 

Z = S + i P = d $£ R , (2.3) 
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wherein we already introduced the angular decomposition of the complex field S into a 
Hermitean field $ and unitary fields ^ and £r which are related to the chiral field by 

U = ^l^R. The latter is conveniently expressed in terms of a chiral angle 

U(x) = exp (iO(x)) . (2.4) 
The quark determinant Af diverges and therefore requires regularization. As in a 



study of the mesonic sector of the NJL model |Tj] as well as in previous studies of the 
soliton sector 0, §, || |10| we will use Schwinger's proper time regularization ]Tj| which 
introduces an 0(4)-invariant cut-off A after continuation to Euclidean space. For the 
regularization procedure it is necessary to consider the real and imaginary part of Af 
separately 

A F = Ar + Aj 

A r = ~Trlog(#£|k) 

Ai = l -Tv\og{{pl)- l p E ). (2.5) 

The real part Ar diverges like logp 2 for large momenta p whereas the imaginary part Ai 
does not contain any divergencies, i.e. it is finite without regularization. However, we 
believe that it has to be regularized also in order to have a consistent model. After all, 
the occurrence of the cutoff is a very crude way of mimicing the asymptotic freedom of 
QCD. 

For the real part the proper time regularization consists in replacing the logarithm by 
a parameter integral 

1 f°° ds / + \ 

Ar — ► —— — Tr exp -sp E pE , (2.6) 



2 Jl/A? s 

which for A — > oo reproduces the logarithm up to an irrelevant constant. Since the 

operator $ e Pe is Hermitean and positive definite this integral is well defined. For the 
imaginary part the regularization procedure is equivalent, 

Ar — -I r -Trexp (s^)- 1 ^) , (2.7) 



2 A/A 2 s 

however, in this case one has to be careful concerning the convergence of the integral, see 
section 3. 

In order to determine the coupling constants g\ and g% from the meson sector it is 
sufficient to only inspect Ar. Varying the regularized effective action with respect to the 
scalar and pseudoscalar fields yields the Dyson-Schwinger or gap equations 

(Sij) = SyM 
M = m°-2g 1 (qq) 

(qq) = -M 3 ^r(-1,M 2 /A 2 ). (2.8) 

The quantity M is the dynamically generated constituent quark mass and (qq) the quark 
condensate. A non-vanishing quark condensate reflects spontaneous breaking of chiral 
symmetry. 



4 



Performing the derivative expansion ||14|| of Ar allows to read off the pion decay constant 



f n as the coefficient of the expression quadratic in the derivatives of the physical pion field 



A R = d 4 x If tr^vr^vr + • • ■. (2.9) 



f 2 

j4„ J_JT_ 

4 

The inclusion of vector and axialvector mesons leads to pseudoscalar-axialvector meson 
mixing, especially tt — at— mixing. This renormalizes the pion field and thus affects the 
pion decay constant which is then given by 



g\ l + QM 2 /m 2 



where 



are the universal vector coupling constant and the vector meson mass. As input quantities 
from experiment we use the pion decay constant f n = 93MeV and the p meson mass 
m p = 770MeV. For a given constituent quark mass M we then determine the cut-off A via 
eqn.( |2.l0| ) and subsequently g<i via fl2.11| ). It is important to note that the 7r — a\— mixing 



increases A significantly, e.g. for M = 350MeV we find A = 1274MeV compared to 
A = 649MeV when ir — a\— mixing is disregarded. Expanding A m up to second order in 
the pseudoscalar fields allows to express the current quark mass m° in terms of g\ and the 
pion mass = 135MeV: m° = gxrn^f 2 /M. Finally we employ the gap equation (|2.8|) to 
eliminate g\ in terms of the constituent quark mass M which from now on is considered 
as the only free parameter of the model. 



3. The Energy Functional for Static Meson Fields 

Next we will consider the energy functional of the static soliton in SU(2). After Wick 
rotation, i.e. xq = —2x4 and Vb = — iVa the Euclidean Dirac operator corresponding to 
eqn. fl2.2|) is given by 

i(3pE = -d T -h, 

h = a-p + iV 4 + ij 5 A 4 + a-V + j 5 a-A + i3{P R E + P L E^) (3.1) 

wherein r denotes the Euclidean time. In Euclidean space r, V4 and A 4 have to be 
considered as Hermitean quantities. This leads to a non-Hermitean Hamiltonian h even 
for static configurations (i.e. [d T ,h] = 0) if non-vanishing time components of vector or 
axialvector meson fields are included. 

We fix the scalar field at its vacuum value $ = Ml. This is a priori an unjustified 
approximation and indeed it has been shown |T9[] that allowing (ft to be space dependent 
leads to a collapse of the soliton in the case when only scalar and pseudoscalar fields are 
present. Employing, however, a well motivated generalization |2(J of the NJL model which 
includes an additional $ 4 potential in the mesonic action ( |2. 2| ) yields a stable soliton with 
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the scalar field deviating only slightly from its vacuum value pl|, p^j . Thus it is reasonable 
to keep the NJL model (|2~2|[2.5| ) and restrict the scalar field to its vacuum value. 
For the chiral field we impose the hedgehog ansatz: 

U{x) =exp(iT-f0(r)). (3.2) 

This configuration has vanishing 'grand spin' G = I + cr/2 + r/2, i.e. [G, U] = 0. The 
only possible ansatz for the isoscalar-vector field to with grand spin zero has vanishing 
spatial components^ = 0): 

V* = w M = w(r) V (3.3) 

Parity invariance requires the isoscalar-axialvector meson field to vanish in the static 
limit. For the isovector- (axial) vector meson fields we use the spherically symmetric 
ansdtze 

VI = o, V? = e akl f k G{r), 

A\ = 0, A\ = r¥ a F(r) + 5 ia H (r) (3.4) 

where the indices a,i and k run from 1 to 3. Note that corresponds to the physical 
p meson while the physical axialvector meson a\ is obtained from A® and terms involv- 
ing <9j(f a 0) after removing the pseudoscalar axialvector mixing. The Euclidean Dirac 
Hamiltonian now reads 

h = ct • p + iuj(r) + M/?(cos6(r) + 275T • rsinO(r)) 

+±(a x f)-rG(r) + ±(«r*)(T.*)F(r) + I(<r • r)ff(r) (3.5) 

which obviously is not Hermitean since w(r) is real giving rise to complex eigenvalues of 
h. The matrix elements of the Hamiltonian ( |3.5| ) are evaluated in the quark spinor basis 
proposed in ref . |^3 1 . These spinors are characterized by their grand spin eigenvalue G and 



their parity transformation properties. Since (|3.5|) commutes with the grand spin operator 



and is parity invariant as well, the Hamiltonian ( |3.5| ) is diagonalized for each grand spin 
and parity eigenvalue separately. The momentum eigenvalues of the basis spinors are 
discretized by putting the system into a finite spherical box of radius D and demanding 
the upper components to vanish at the boundary for spinors with parity eigenvalue (— 1) G 
and the lower components for spinors with parity eigenvalue (— 1) G ~ 1 . We list all matrix 



elements of operators appearing in (|3.5|) as well as the explicit form of the basis spinors 



in appendix A. The boundary conditions of ref.|23[ have the advantage that there are no 
spurious contributions to the equations of motion for the u field stemming from the fact 
that we only consider a finite basis in momentum space. The spurious contributions to 



the profile function G(r) of the p field appearing in the basis of ref. |23[ are well under 
control and may explicitly be eliminated. However, we will see later that there are finite 
size effects for the u meson profile originating from a large but finite D. 

For static configurations the eigenvalues of d T , zO n = i(2n + l)vr/T, with {n = 0, ±1, 
±2,..) may be separated rendering the temporal part of the trace feasibleB Thus the 
eigenvalues A n)1/ of the operator d T + h read: 



K,u = -i£l n + t v = -iQ n + + ie J u . (3.6) 



*The eigenfunctions of d T assume anti-periodic boundary conditions in the Euclidean time interval T. 
The Q n are the analogues of the Matsubara frequencies with T figuring as inverse temperature. 
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The fermion determinant is expressed in terms of the eigenvalues X n y. 

i fi =Vlog(A„,A;j and A T = Wlog^) . (3.7) 

u,n u,n ^n,v 

Using ( |3.6| ) the real part reads: 

A R = ^£log((O n - e[f + (a 2 ) 



according to the proper time regularization scheme (|2.6| ). For large Euclidean time inter- 
vals (T — > oo) the temporal part of the trace may be performed 

where we have shifted the integration variable z — el — > z. For T —>■ oo we may read off 
the Dirac sea contribution to the real part of the energy functional from Ar — > —TE^ ac : 

^c = T%El^ir(-^(^/A) 2 )- (3-10) 

V 71 " v 2 

For the imaginary part we obtain 

^ = ~(£ £ i°g(A,, B )-£ £ MA:,j) = iE £ lo s|r3 ( 3 - n ) 

Z k n=-oo v n=-oo Z f n=-oo 41 4 " e i/ 

where we have reversed the sign in the first sum over the integer variable n. Next we 
express Ai in terms of a parameter integral: 

1 00 r 1 —if 1 
^= EE /> .-n ^77 - ( 3 - 12 ) 



In analogy to ( ft.lOQ we may carry out the temporal trace in the limit T — ► oo: 

^=T?/> r £^ [<( ,-4-^ - (313) 

Shifting the integration variable z — \t{, —* z the integral over A may be done 

* = T?^££?W (3 - 14) 

^4.7 is regularized in proper time by expressing the integrand as a parameter integral: 

\ roc 

^.^tepl-rl^fffl (3.15) 



which obviously is finite for A — > oo. Continuing the evaluation of Aj in analogy to 
eqs. (|3.8[ - |3.ld|) we find for the contribution of the Dirac sea to the imaginary part of the 
Euclidean energy E^ ac 

^ ac = Z ^E4sign( e f){ 1 A > A / not ^larized (3.16) 

vac 2 y lAv, Aj regularized v ' 

where 

N v = -^=r(^, (e^/A) 2 ) (3.17) 

are the vacuum "occupation numbers" in the proper time regularization scheme. The 
upper case, of course, corresponds to the limit A — > oo. Obviously only the real part of 
the one-particle energy eigenvalue is relevant for the regularization of Ai- Eqn. ( |3.16| ) 
reveals that we have succeeded in finding a regularization scheme for Aj that only involves 
quantities which are strictly positive definite. This is not evident from the definition of 

Aj(B- 

For soliton configurations with vanishing lj (i.e. e u = e„) there is no contribution from 
the imaginary part and eqn. (|3.10|) is the expression for the energy of the Dirac sea. 
The total energy functional contains besides E^ c and E^ &c also the valence quark energy 

*W = N c J2Vv\e?\ E 1 ^ = N c J2v»sign(e?)el (3.18) 

V V 

with r)n = 0, 1 being the occupation numbers of the valence quark and anti-quark states. 
Furthermore the meson energy is obtained by substituting the ansatze (|3.2| - |3.4j ) into ( |2.2|) : 



E m = 4vr J drr 2 (ml / 2 (1 -cos 9(r)) 



+ (!^) 2 ( G 2 (r) + l F 2 (r) + F(r)if ( r) + 3^ 2(r) _ 2u j 2 (r))) (3.19) 
g v 2 2 ' 

Note that we are working in the isospin limit which implies = m p . Continuing back 
to Minkowski space we find for the total energy functional: 

E[Q, u , G, F, H] = + E^ + E* c + E^ c + E m . (3.20) 

The equations of motion for the meson profiles are obtained by extremizing the static 
Minkowski energy (|3.20|) . In a generic way we may write: 

wherein denotes any of the meson profiles 0, G, u, F or H . Since h is not Hermitean 
(in Euclidean space) we have to distinguish between left and right eigenvectors of h. The 
corresponding eigenvalue equations read: 

fr|tt„) = e„|tt„) l$t v \h = e v ($! v \ i.e.tf\V v ) = (Z\* v ). (3.22) 

The normalization condition is (^f u ) = 5^ u . In order to evaluate the derivatives 
8e*/5(j) it is helpful to decompose the Hamiltonian operator ( ^.5[ ) into Hermitean and 
anti-Hermitean parts 

h = h e + iu (3.23) 
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where h® includes all Hermitean terms of the Euclidean Dirac Hamiltonian ( |3.5| ). Obvi- 
ously both, h® and u, are Hermitean implying \^ u ) = We may therefore extract 
the real and imaginary parts of the one particle energy eigenvalue 

^ = \{(K\h\*v) + 

4 = liiKWJ - {*v\tf\K)) 

= (*?\u\V?) - (*l\u\*l) + (*l\h e \yX) + (^\h e \^ u ) (3.24) 

where we employed the decomposition \ty u ) = |\E^) + Note also that = 

(^^1 + We are now equipped with expressions for the real and imaginary parts of 

the energy eigenvalues e M which are suitable to evaluate the derivatives with respect to 
the meson fields. For example we have: 

— r 



8u(r) 



/^((r|^)(^|r)-(r|^)(^|r)). (3.25) 



The expressions for the functional dependence of the energy eigenvalues may now be 
substituted into the equations of motion ( |3.21| ). The individual equations for Q,u,G,F 
and H are displayed in appendix B where we also discuss the spurious contributions to 
the equation of motion for G. 

At this point it is indispensable to explain the differences between our approach and 
the treatment in ref. [piR . We would like to stress that in order to regularize the fermion 
determinant a continuation to Euclidean space has to be performed. In Euclidean space 
the regularized fermion determinant is well defined in terms of the complex energy eigen- 
values e M of the Euclidean Dirac Hamiltonian ( p . 5|) . This final expression for the fermion 
determinant may then be continued back to Minkowski space yielding the energy func- 
tional( F20|) - Due to the non-linear structure of the regulator functions in ( |3.10|) and 



(|3.16| ) the operations "regularization" and "continuation" do not commute. However, 



this is just the procedure used by the authors of ref. p4|| : The energy eigenvalues had 
firstly been continued to Minkowski space and then the energy functional was regularized. 
In appendix C we present a detailed comparison of the two different treatments with the 
help of a toy model. 



4. Numerical Results 

In this section we present the numerical results characterizing the soliton solution. 
First we will explore the soliton containing all meson fields as discussed in the preceding 
sections. These results will then be compared to cases where some of the vector meson 
fields are switched off in order to examine the effects of various vector meson fields on the 
soliton. 

The meson profiles of the soliton solution are determined by iteration. /. e. we start off 
with test profiles for 0, u, G, F and H to evaluate the energy eigenvalues e M of the static 
Hamiltonian ( |3.1| ). Subsequently the meson profiles undergo modifications according to 
the equations of motion ( ft.21| ). In turn these modified profiles serve as input for the 
static Hamiltonian. This process is repeated until a convergent solution is obtained, i.e. 
a self-consistent meson field configuration is constructed. 
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Table 4.1. The soliton energy E as well as its Dirac sea and mesonic contributions E vac and E m for 
different values of the constituent quark mass M . For the pion mass 771 ^ the physical value (135MeV) 
and the chiral limit (777^ = 0) are considered. Also shown is the energy of the 'dived' level (e va i). 



M (MeV) 


300 


350 


400 


777,,- (MeV) 





135 





135 





135 


E (MeV) 


1117 


1155 


1024 


1061 


952 


980 


EL (MeV) 


690 


690 


569 


572 


504 


503 


^vac (MeV) 


36 


36 


31 


30 


24 


24 


E m (MeV) 


391 


429 


424 


459 


423 


453 




-0.11 


-0.11 


-0.44 


-0.42 


-0.60 


-0.61 


<JM 


0.10 


0.10 


0.09 


0.09 


0.07 


0.08 



Table 4.2. Same as table (4.1) for 777^ = 0, however the imaginary part is not regularized. 



M (MeV) 


300 


350 


400 


E (McV) 


1227 


1210 


1175 


EL (MeV) 


669 


563 


460 


EL (McV) 


55 


257 


298 


E m (MeV) 


335 


389 


417 


<&JM 


0.02 


-0.28 


-0.47 




0.13 


0.16 


0.16 



We have found stable self-consistent solutions for constituent quark masses in the range 
300MeV< M <4OOMeV0 In all these calculations the parameters have been kept at their 
physical values (cf. section 2). The existence of these stable solutions represents the main 



difference as compared to the results of ref.[24]. Those authors only find solutions for 
m w > 870MeV. 

In table (4.1) we display the energy E of the self-consistent soliton solution. This table 
also contains the various contributors to E as they appear in eqn.( |3.20 ) for several values 



of the constituent mass M. Furthermore we compare to the results obtained in the chiral 
limit (m^ = 0). 

The most striking result observed from table (4.1) is the fact that the real part of the 
energy eigenvalue associated with the valence quark state is negative! I.e. the valence 
quark has joined the Dirac sea and thus the baryon number is completely carried by the 
polarization of the Dirac sea. Thus the soliton of the NJL model supports the picture of 
baryons as topological solitons of meson fields. As a reminder we would like to mention 
that for eL < the valence quarks' contribution to the energy is already contained in EL 
and EL- Thus it must not explicitly be added in ( |3.18| ). We also observe from table (4.1) 
that if the imaginary part of the fermion determinant is regularized it only contributes a 
minor part to the total energy. However, we see from table (4.2) that this is no longer the 
case when the regularization of the imaginary part is abandoned. In that case the total 
energy E is almost independent of the constituent quark mass M while E decreases by 
about 20% in case the imaginary part is regularized when M is changed from 300MeV to 
400MeV. For a non-regularized imaginary part the valence quarks appear to be slightly 



*We do not exclude the existence of stable solutions for an even larger range of M . 
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Table 4.3. The soliton energy for various treatments of the NJL soliton. The meson fields in the first line 
denote the allowed meson profiles. All numbers are evaluated for a constituent quark mass M = 350MeV 
and = 0. If present, the imaginary part is regularized. 





IT 


71, p 


7T, UJ 


ir,p, ai 


ir,UJ,p, ai 


E (MeV) 


1214 


957 


1310 


1010 


1024 


EL (MeV) 


561 


655 


629 


615 


569 


EL (MeV) 








-56 





31 


E m (MeV) 





155 


-42 


395 


424 


Ci/M 


0.62 


0.14 


0.50 


-0.13 


-0.44 


<JM 








0.24 





0.09 



above the Dirac sea for M=300MeV. For all other cases the valence quarks' energy remains 
negative, although less strongly bound in the Dirac sea than for a regularized imaginary 
part where the uj profile is less pronounced. Table (4.1) also reveals that the inclusion of 
a finite pion mass does not alter the results for the energy drastically. The total energy 
is increased by about 40MeV due to the additional term in E m fl3.19| ). 

Next we wish to investigate the role of the different meson fields for the NJL soliton! 
In order to do so we compare in table (4.3) the results obtained for the soliton energy in 
cases with different (axial-) vector mesons incorporated. Obviously, the inclusion of the 
uj field always increases the soliton energy while the p and a± lower the energy. Though 
the latter result is anticipated the former is somewhat surprising since the meson profiles 
are obtained by extremizing the total energy. However, this increase is understood by 
taking into account that the uj field is proportional to the baryon number density which 
in turn is constrained by unit baryon number. We also observe from table (4.3) that the 
inclusion of any of the (axial-) vector mesons lowers the energy eigenvalue of the valence 
quark. The a\ obviously effects the valence quark to join the Dirac sea. This import 
result was already obtained previously [|TT|] and the main conclusion to be drawn out of 
our present calculation is the fact that this result is not spoiled by the to meson which 
was not included in ref.|TT[; on the contrary the uj meson gives an even stronger binding 



of the valence quark. This may be understood by noting that the uj is repulsive yielding 
a large spatial extension of the soliton. This, in due, causes the valence quark energy to 
drop. 

The repulsive character of the uj field may also be observed directly from radial behavior 
of the chiral angle, G(r). In fig. (4.1) we display 0(r) for various treatments of the NJL 
soliton. G(r) develops the largest tail in the case when the uj meson field is the only 



one added to the chiral field|12|. The axialvector field provides a significant attraction 
resulting in a slope for the chiral field which is larger than in the case when G(r) is the 
only field being present. The inclusion of the uj meson on top of the isovector mesons p 
and ai alters the chiral angle only slightly. A similar behavior can be observed for the p 
meson profile G(r) (cf. fig. (4.2)) as well as the axialvector meson profiles F(r) and H(r) 
(cf. fig. (4.3)). On the other hand the inclusion of the isovector mesons on top of the 
7i — uj system significantly reduces the strength of the uj meson profile as may be seen in 
fig. (4.4). 

Of course, in the present form the soliton carries neither good spin nor isospin quantum 
numbers, hence, it does not describe physical baryon states. Following the work of Adkins, 



^We do, however, not consider a system containing the p and uj mesons besides the pseudoscalar fields. 
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Fig. 4.1. The chiral angle 0(r) as a function of the radial distance T. The solid line corresponds to the 
case when all vector mesons are included; the dashed line to the 7T — UJ system; the short dashed to the 
7T — p system; the long dashed to the 7T — p — d\ system and the dashed dotted to the case when the 
chiral angle is the only field present. 



Fig. 4.2. The vector meson profile G(r) as a function of the radial distance r. Solid line: all vector 
meson fields are present; short dashed denotes the 77 — p system and the long dashed to the 7T — p — d\ 
system. 



Nappi and Witten|25|] the soliton is projected onto baryon quantum numbers by a cranking 
procedure which introduces time dependent collective coordinates for the zero modes R(t). 
Explicitly one imposes the ansatz: 



£(x,t) = R(t)£(x)ltf(t) R(t) G SU(2). 



(4.1) 



and similarly for the isovector fields p and a\. Although this approach has been suc- 
cessfully applied to the SU(2) NJL model of pseudoscalar fields only|26| and even been 
generalized to SU(3)[|27j this is not the whole story in the presence of vector mesons. 
The collective rotation R(t) excites additional vector meson components as e.g. the space 
components of uo and the time components of p which are absent in the static caseQ. 
Since the investigation of these excitation is beyond the scope of this paper we will not 
discuss the resulting baryon spectrum, especially the nucleon-A mass splitting. Neverthe- 
less there are a few baryon properties which dependent only on the static fields and may 
thus be evaluated at the present stage; the most prominent of which is the axial charge 
of the nucleon, qa- 



In the NJL model the current field identities hold|J14|. Therefore the axial current J, 



is directly proportional to the axial field: 



5/j 



J< 



5/j 



(4.2) 



wherein the superscript denotes the isospin component. Noting that gA is obtained as the 
matrix element of at zero momentum transfer we immediately obtain 



9a = — 



2tt 

92 



drr z 



H(r) + -F(r) 



(D 33 ) p . 



(4.3) 



{Dszip = 1/3 refers to the matrix element of ti(Rr 3 RW 3 ) between proton states. It is 
important to mention that eqn. (|4.3|) represents an exact result which is not subject to 



Fig. 4.3. The axialvector meson profile F(r) and H[r) as functions of the radial distance r. H[r) 
is non-vanishing at the origin. The case when all vector mesons are present is denote by the solid and 
dashed lines. The 7r — p — CL\ system is represented by the long dashed and dashed dotted lines. 
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Fig. 4.4. The vector meson profile U)(r) as a function of the radial distance r. The solid line denotes the 
case when all vector mesons are present while the dashed line represents the 71 — UJ system. 

Table 4.4. The axial charge Qa of the nucleon in the various treatments of the soliton in the NJL model. 
The last line corresponds to the case when the imaginary part is not regularized. 



M (MeV) 


300 


350 


400 


71 




0.78 


0.73 


7T, UJ 




0.98 


1.03 


n,p,ai 


0.31 


0.27 


0.13 


n,uj,p,ai 


0.34 


0.25 


0.23 


7r,u,p,at 


0.54 


0.39 


0.28 



renormalization due to tt — a\ mixing. Making use of the equation of motion for the 
axialvector profiles H(r) and F(r) (see appendix B) we may reexpress gA as a mode sum 
over quark spinors: 

9A = -^E{ [«ksT3|V?> - WksTsIt© + «k3T 3 |V>;> +<V>>3T 3 |^)b, 

+ Wk3T3|V?)]//(WA) } (4.4) 

The regulator functions fnj are listed in appendix B (cf. eqns. ( B.3| , |B.4| )). The mode sum 
( |4.4| ) has the advantage that it may be employed in models without axialvector mesons. 

Unfortunately our numerical results which are listed in table (4.4) are somewhat dis- 
couraging since they are well below the numerical value = 1.25. This is especially 
pronounced in case the valence quark energy is negative. In order to understand the 
origin of this shortcoming let us consider the case when the isoscalar field uj is absent. 
Then eqn. ( |4.4j ) reduces to 



N c N c 
9a = — ^-^vai(val|o- 3 r 3 |val) + — ^ sign(e / ,)erfc(|e At |/A)(/x|cr 3 r 3 | y u) (4.5) 



which is identical to the expression derived in ref.[2S]. Eqn. ( |4.5| ) reveals that once the 
valence quark energy has become negative its contribution to gA is strongly suppressed 
by the proper time regularization. This suggests that a regularization prescription which 
does not affect low-lying states as strongly as the proper time scheme would be highly 
desirable in order to describe correctly. This consideration is supported by a simple 
modification of eqn. (|4.5|) . In the sum over all eigenstates we replace the complementary 
error function by a sharp cut-off function. Then the contributions from the low-lying 
states are not affected by the regularization procedure. Choosing, e.g. a constituent 
mass of 300MeV the prediction for increases drastically to 1.04 in the tt — p — ai 
system. Of course, this exploration does not represent a consistent calculation but merely 
demonstrates that strongly depends on the regularization description. If we consider 
the 7i — uj system the strong repulsion of the uj field transfers to an increased prediction 
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Table 4.5. The isoscalar radius (r 2 = ^) l l 2 of the nucleoli in various treatments of the NJL model. The 
radii are given in fm. 



M (MeV) 


300 


350 


400 


7T 




0.89 


0.76 


7T,p,ai 


0.55 


0.55 


0.56 


TV, UJ 




2.06 


1.77 


7r,u,p,at 


1.39 


1.40 


1.29 



for gA- This is also obvious from fig. (4.1) since in the absense of axialvector fields qa may 



also be obtained from the size of the "pion tail" ^5 . 

We may also investigate the isoscalar radius of the nucleon without explicitly perform- 
ing the collective quantization. Again due to the current field identity the isoscalar charge 
density is proportional to the uj profile yielding the isoscalar radius: 

4m 2 

N^gJ 

The normalization factor 



4?7l 2 

(r 2 =0 ) = N" 1 — I d 3 r r 2 u{r). (4.6) 



Am 2 r 

is one if the imaginary part of the determinant remains unregularized. 

As already indicated in the previous section our numerical results contain some spurious 
finite size contributions. These finite size effects are manifested in a contribution to the uj 
meson profile which is proportional to D~ 3 . Especially we find that our numerical solution 
for u(r) accquires a finite value at the edge of the box: u(r — D) ~ D" 3 . From eqn. 
(|4.6j ) it may be observed that then (r] =0 ) would diverge as D — > oo. We eliminate this 
spurious contribution by enforcing u(r) to vanish at the boundary. This is accomplished 
by including a Lagrange multiplier A in the energy functional ( [3.20| ): 

E[Q, uj, G, F, H] — ► E[Q, uj, G,F,H} + \ J d 3 rf e (r)uj 2 (r) (4.8) 

wherein f e (r) is a positive radial function which vanishes everywhere except within a 
small vicinity e of r = D. By iteration of the modified equations of motions A is adjusted 
such that the additional term vanishes. Of course, there is some arbitrariness in choosing 
/ e (r). We demand the change of the total energy with A included to deviate from the 
case A = by less than lMeV. In the region of physical interest (r < 2fm) u(r) remains 
almost unaltered and we obtain the desired effect that (r| =0 ) stays finite as D — > oo and 
assumes a constant value. In case the uj field is not present we use the corresponding 
unregularized mode sum ( |B.6| ) as input in eqn. Q4.6|). The resulting data are displayed 
in table (4.5). The repulsive effect of the isoscalar vector field uj is obvious and is not 
completely compensated by the attraction provided by the isovector fields p and a±. I.e. 
the prediction still overestimates the experimental value (r| =0 ) 1//2 w 0.8fm. 



5. Conclusions 
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We have found stable soliton solutions in the NJL model with all low-lying vector and 
axialvector meson fields included. The isoscalar vector meson u is excited only if the 
imaginary part of the Euclidean action is taken into account. It is essential to note that 
the necessary regularization can only be performed in Euclidean space. This implies that 
the energy functional can be defined properly only if one continues to Euclidean space, 
regularizes and then continues back to Minkowski space. This is quite distinct from the 
procedure of ref. |2l| where the one-particle energy eigenvalues had firstly been continued 
to Minkowski space and the resulting energy functional was "regularized" . As a matter 
of fact, their results are quite different from ours. Once they incorporate the u meson 
they do not even find stable soliton solutions for the physical value of the us mass. 

For reasonable values of the constituent quark mass we find that the energy eigenvalue 
corresponding to the valence quark state is negative. Therefore the baryon number is car- 
ried by the asymmetry of the vacuum. Thus the NJL model supports Witten's conjecture, 
i.e. the Skyrmion picture of the baryon. In this respect we remark that the u meson is 
even decreasing the valence quark level despite its repulsive nature. This can be under- 
stood be noting that the u meson broadens the other meson profile functions thereby 
increasing the binding energy for the valence quark state. However, if the imaginary part 
is not regularized this effect is inverted. 

The calculation of the isoscalar mean square radius of the nucleon is plagued by a 
numerical deficiency stemming form the finite spatial extension where we iterate the 
equations of motion. We have introduced a Lagrange multiplier to keep the result for 
the isoscalar radius finite. The prediction for the isoscalar mean square radius of the 
nucleon overestimates the experimental value by about 40-50%. Nevertheless we have 
seen the repulsion provided by the isoscalar vector field as well as the attraction due to 
the isovector (axial-) vector fields. 

The axial coupling qa comes out much too small. We have demonstrated that this 
behavior is due to a deficiency of the proper time regularization which attaches "weight 
factors" smaller than unity already to levels with quite a small energy. We therefore 
conclude that a regularization which does not effect low energy levels is highly desirable. 

In order to calculate other observables one has to project the soliton solution on good 
spin and flavor states as was done for the NJL soliton with pseudoscalars only in refs. |26|, 



29], |3(|. However, including the vector and axialvector mesons this is a quite involved 
task since components which vanish on the classical level but get excited by the collective 
rotation (4.1) have to be taken into account. On the other hand, we have seen that the 
NJL soliton is close to a Skyrmion. Thus it makes sense to use the self-consistent meson 
profiles in an action obtained by the gradient expansion of the fermion determinant. A 
first step towards this direction is to calculate static properties in this approximation and 
compare those with the exact results. 

Finally we would like to remark that problems arising from the regularization can be 
overcome only if one takes a renormalizable quark model as starting point. E.g. one 



way were to use a bilocal effective quark interaction leading to bilocal meson fields|31 



Whether such models allow for soliton solutions and whether these solutions support 
Witten's conjecture is, of course, an interesting question which deserves further studies. 



Appendix A: Matrix elements of the static Hamiltonian 

For our numerical calculations we have used as basis the orthonormal eigenstates of 
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the free Dirac Hamiltonian 



ho = ap + f3M. 



(A.l) 



Since the static Hamiltonian commutes with the grand spin operator G = I + cr/2 + r/2, 
our basis spinors are characterized by the grand spin quantum number G in addition to 
the total angular momentum j = G±l/2 and the orbital angular momentum / = G, Gil. 
The latter also determines the parity eigenvalue of the spinor, tt = (—1)'. The grand spin 
eigenstates: 



\GMjl) 



m,j ,m,[ 



{ ~ y jm J ,l/2M-mj i ~'lm l ,l/2m : j-mi 



1 1 

\lmi)\-mj - mi) s\-M - rrij) i (A. 2) 



are two component spinors in both spin(S')- and isospin(/)-space. The momentum, and 
therefore the energy, is discretized by putting the system in a spherical box of radius 
D and requiring appropriate boundary conditions. For the calculations reported in this 
paper we used the boundary conditions of ref. P^fP The discrete momenta are then 
given by requiring 



j a (<gD) = 0. 



(A.3) 



where ja are the spherical Bessel functions. The energy eigenvalues corresponding to 

the momenta q n in (|A.3|) are = ±\J (<?^) 2 + rn? . Now the Dirac spinors are easily 
constructed: 



(i) J = G + \- l = G: 



\G + + n) 



AC 



G 

G+ln 



^< + J G (q^r) \G M (G + |) G) 
wC-j G+ i(q%r) \GM{G+\) (G+l)> 



(A.4) 



(ii)j = G-±; l = G 



lr _ , , KfG ( iwC + j G (qC r ) \GM{G-\)G) 

' } " G+ln [ -<-JG-i(q°r) |GM(G-§)(G-1)> 



(A.5) 



(hi) J = G+\- l = G + l: 



\G + -n) 



G+ln 



\G M (G+|) (G+l)) 
\G M (G+§) G) 



(A.6) 



(iv) j = G-h l — G—1: 



2' 



\G 



n) 



G+ln 



iw^ + j G ^r) \GM{G-\) (G-l)) 
vfiTMr) \GM{G-\)G) 



(A.7) 



*For a discussion of a different type of boundary condition see ref. 
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where we used the following abbreviations 



Gn 
L+ 



D-S\3a{&D)\-i 



'1 + 



m 



rn 



(A.f 



We start with listing a few helpful relations involving the two-component grand spin 
eigenstates ( |A.2| ) only: 



(GMhhlT-rlGMM 



i 



2G+1 

-1 
2G+1 



for ji = j 2 = G + ~ , \l 2 



for ji = j 2 = G 



2 ' 



\k-k\ 



2 V /G(G+1) 
2G+1 





for h = G±i 



2' 



J2 



Gt«, I ^2 — ^1 1 — 1 



(A.9) 



otherwise 



(GMfrhKtr xf)r\GMj 2 l 2 ) 



* + + (A.10) 

-J2O2 + 1) + k(k + MGMjMt ■ i\GMj 2 l 2 ) 



{GMjih\(cr-r)(T ■ i)\GMj 2 l 2 ) 



and 



-1 



2G+1 

1 

2G+1 

2 V /G(G+1) 
2G+1 



2 V /G(G+1) 
2G+1 



for 
for 

for 
for 

otherwise 



Ji 
Ji 

h 
h 



j 2 — G + 2 , h — k 
ji — G — I ,h=k 



Gi;, j 2 — G ^ 5 



2 ' 

G 



2' 



(All) 



Ji 
\l2-h 



G±h 32 



2 ■ 
2 



Gt\ 



2' 



1 


for ji 


= 32 = 


G± 2 


k = 


k 


= G±l 


2G-1 
2G+1 


for ji 


= 32 = 


G - 5 


h = 


k 


= G 


2G-3 
2G+1 


for ji 


= 32 = 


G + 5 


h = 


k 


= G 


4 V /G(G+1) 
2G+1 


for ji 


= G± 


i 32 = 


Gt 


1 

2' 


h = k = 



(GM 3l h\(a ■ r)\GMj 2 l 2 ) 



otherwise 
A few helpful identities to verify the above expressions are: 



(A.12) 



G 



[(o-f)(r-f)f 
[<x • r - (<x-r)(T -r)f 



4 
1 



2-2(o--f)(r-f). 



(A.13) 
(A.14) 
(A.15) 
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In the following we shall display all matrix elements which are needed for the static 
Hamiltonian. 

We introduce some additional abbreviations: 



1. 6 = Pf(r 



nm 

IV 

nm 
IV 



N„ 



Wn~w^-ji{q^r)j v (qgr) 
Wn~wg + ji{q%r)j v (qgr) 
\3G + i(qZD)nj G+1 (qgD)\-i 



(A.16) 



(G + + n | O \G + + m) 
(G- + n | 6 \G- + m) 
(G + - n\6 \G + - m) 
(G n\6 \G m) 



N nm Jo dxX 2 f(Dx) (a>nm ^nm 1G+1 

N nm Jo dxX 2 f(Dx) (dnm ~~ ^nm^ 

N nm Jldxx 2 f{Dx) (a%+ 1G+1 

N nm Jldxx 2 f(Dx) (a " 10 - 1 



nm ) 
b GG ) 

nm w nm J 



2.0 = it- f/3 75 /(r) 



(G + + n 


6 


\G + 


+ m) 


(G- + n 


6 


\G- 


+ m) 


(G + + n 


6 


\G- 


+ m) 


(G- + n 


6 


\G + 


+ m) 


(G + -n 


6 


\G + 


— m) 


(G--n 


6 


\G- 


— m) 


(G + -n 


6 


\G- 


— m) 


(G--n 


6 


\G + 


— m) 



M+1G 



3. 6 = (a x f)-T/(r) 



N nm Jq dxX f(DXj of;-!-.!. 1 ' nm 

N nm J 1 dxx 2 f(Dx^ - ( r-G—}< 



^mn I 



2G+1 y C nm 



+ C 



G-1G\ 



j dxx 2 f(Dx) — 2G+1 — 



M-lG 



c 



G+1G\ 



N nm J 1 dxx 2 f(Dx) 



2G+1 

-1 



GG+1 i GG+1\ 

™ "T 771.71 J 



2G+1 v^nm 



- 1 (c^- 1 + c^- 1 : 



^nm/o d.r.r ./ i /<*•''/ 2(V+! " ""' 
N nm Jo 1 dxX 2 f(Dx) VG(g+l) ^ GG -! _ ^GG+1 

N„ m Jo 1 dxx 2 f(Dx) 



2G+1 v^nm ''mn 



2G+1 



„GG-1 _ „GG+1 
^mn ^nm 



^We employ the standard Dirac representation for the 7 matrices. 
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(C + + n 
(G- + n 

(G + + n 

(G-+n 
(G + -n 
(G--n 

(G + -n 

(G--n 

4. 6 = (<r-f)(r 

(G + + n 
(G-+n 

(C + + n 

(G- + n 
(G+-n 
(G--n 

(G + -n 

(G--n 



O \G + + m 
6\G-+m 

6 \G-+ m 

6 \G + + m 
6\G + -m 
6\G--m 

6\G--m 

6 \G + - m 

f)/(r) 

6 \G+ + m 
6\G- + m 

6\G- + m 

6 \G + + m 
6 \G+ - m 
6\G--m 

6\G--m 

6\G + -m 



N nm jQ<ixx 2 f{Dx) 
N nrn ^dxx 2 f{Dx) 

N nm J l dxx 2 f(Dx) 

N nm fo dxX 2 f(Dx) 

N nm ^dxx 2 f{Dx) 
N nm j l dxx 2 f(Dx) 

N nm fo dxX 2 f(Dx) 

N nm J l dxx 2 f(Dx) 



2(G+1) / G+1G I r .G+lG) 

2G+1 \ nrn ' ^mn ) 



-2G ( G-1G I G-1G\ 



2G+1 ^ nm 



2VG(G+1) ( G-1G _ r G+lG\ 
2G+1 v ran °nm I 



2 V /G(G+1) / G-1G _ r G+lG 
2G+1 \ ram mra 

) 



) 



2(G+1) / GG+1 , „GG+1 
2G+1 \ nrn ~ ^mn 



-2G (JSG-1 



GG-l 



) 



2 V /G(G+1) , GG+1 _ „GG-1 



2G+1 



2 V /G(G+1) , GG-l _ GG-1\ 
2G+1 V nm °mn J 



N nm ^dxx 2 f{Dx) 
N nm J^dxx 2 f(Dx) 

N nm Jo dxx 2 f(Dx) 

N nm j^dxx 2 f(Dx) 
N nm J^dxx 2 f(Dx) 
N nm ^dxx 2 f{Dx) 

N nm J l dxx 2 f(Dx) 

Nnmlo dxX 2 f(Dx) 



_zJ_f a GG, uG+lG+l) 
2G+1 V nm ' u nm ) 

-J—(a GG + hG-lG-l) 
2G+1 V nm ~ u nm ) 

'a 



2y/G(G+l) ( n GG _ lG+1G-1\ 
u nm I 



2G+1 



2 VG(G+1) , QQ _ h Q_ 1G+1 s 

2G+1 \ u nm u nm ) 

+ b GG ) 
1 u nm I 



-1 I G+1G+1 
2G+1 V nm 

1 <„G-\G-\ 



2G+1 V a nm 



+ 



2 V G(G+1) / Q+iG-i 
2G+1 V u nm 



2 V /G(G+1) / G-iq + i 



2G+1 



- e) 



5. 6 = (<r • r)/(r) 

(G + + n | 6 |G + + m) 
{G - + n\6 \G - + m) 

{G + + n\6 \G- + m) 

{G- + n\6 \G + + m) 



= N nm f Uxx 2 f(Dx)(- 2 ^a GG + b G +^) 



N nm ^dxx 2 f{Dx) (- 



2G+1 

2G-1„GG 



2G+1 a nm 



+ foG-lG-n 



N^dx^Dx) (±*gpO 



N nm J^dxx 2 f(Dx) ( 



2G+1 

4 V /G(G+1) 
2G+1 ' 



,GG\ 



19 





n 


6 




m) 


= N nm jldxx 2 f{Dx) 


( G+1G+1 _ 


2G+3 lGG 
2G+1 nm 


(G - 


n 





\G 


m) 


= N nm ^dxx 2 f(Dx) 




2G+1 ™™ 


(G + - 


n 


6 


\G- - 


m) 


= N nm jldxx 2 f(Dx) 


, Vg(G+1) 
v 2G+1 


nm / 


(G-- 


n 


6 


\G + - 


m) 


= N nm J^dxx 2 f{Dx) 


, Ay/G{G+1) 
v 2G+1 


^nm ) 



Appendix B: Equations of motion for the meson profiles 

In this appendix the equations of motion (|3.21|) are displayed in detail. The functional 
derivative of the one particle energies with respect to the fields yields expressions involving 
the eigenf unctions ( |3.22|) . On the other hand, the derivative of the total energy with 
respect to the one particle energies results in regularization functions for the vacuum part 
and the occupation number for the valence part. As the eigenfunctions (|3.22|) occur only 
in certain combinations in the equations of motion it is convenient to define quark density 
matrices. 

The quark scalar density matrix p(x, y) may be decomposed into contributions due to 
valence and sea quarks: 



p(x,y) 


= p^ + pT + pT 


+ pT 












= E{«(» 






y) 


}Vu, 




P7 l (*,y) 


V 






y) 


}Vu, 






V 


)-ti 




[v] 


}fR 


%/A) 


P? c (*,y) 








y) 


}//( 


e v /A). 



(B.l) 

Similarly, the quark number (or baryon number) density matrix b(x, y) is written as 
b(x,y) = bf + bf + bT + br, 
b™\x,y) = E(»(»)-«(»)K. 

V 

bT\x,y) = -J2{^(x)^(y)+^(x)^J(y)} Vu , 

V 

bT c (x,y) = -J2{^(x)^(y)+^(x)^ I J(y)}f R (e 1/ /A). (B.2) 

V 

Note that the baryon density matrix b(x, y) differs from the scalar density matrix p(x, y) 
not only by the additional factor 7 but also by an exchange of the regulator functions 
fji and fj which are given as the derivatives of the energy functional with respect to the 
energy eigenvalue e u : 

|sign(e^)AC Ai not regularized, 

f R (e u /A) = { ' (B.3) 

-isign(e^)AC + ^(e^/A)exp(-(e^/A) 2 ), Ai regularized 
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fi(eu/A) 



-sign(e^ 



1, Ai not regularized 
Af v Ai regularized 



(B.4) 



with the vacuum occupation number M v being defined in ( |3.17| ). With these definitions 
the equations of motion (|3.21 ) for the meson profiles G, lu, G, F and H read 



sin 6 (r 
oj{r 
G{r 
F(r 
H(r 



M 



Ami J 47T 



iV c tr J — (sin O(r) — ij 5 T ■ f cos Q(r)^jp(x, x), 

r JO 

b(x,x), 

9v N c ti y^((7Xf).T)p(aj,aj), 



Am 2 , 



9v 



N c tr J ^ pfa ■ f)(r ■ f) - (a ■ r))p(x, x) 



Am 2 p 



9v_ 

4m^ 



A^ c tr J ^ P[(cr ■ r)(r ■ r) ~ (a ■ T ))p(x,x). 



(B.5) 
(B.6) 
(B.7) 
(B.8) 
(B.9) 



The traces are over Dirac and isospin indices only. Note that only the u meson profile 
is given by a trace over the baryon number density. All other meson profiles have as 
"source" the scalar quark density. Without regularization the difference between both 
densities would have only been a factor 70 = (3. However, the regularization makes the 
relation between the two densities, and therefore the relation between the u and the other 
meson profiles, highly non-linear. 

For the basis of ref . |23| the RHS of eqn. ( p.7|) does not vanish for the vacuum config- 
uration yielding some spurious contributions to G(r). These may, however, be eliminated 
explicitly by subtracting the corresponding value of the RHS of eqn. ( |B.7| ) at each step 
of the iteration procedure. These spurious contributions are absent when the boundary 
conditions proposed in ref. |27[ are employed at the expense of a non-vanishing RHS of 
eqn. ( |B.6| ) for the vacuum configuration. This then would give similar spurious contri- 
butions to u{r). Thus we conclude that for the problem at hand both sets of boundary 
conditions are equally well suitedu 



Appendix C: Continuation to Euclidean space in a toy model 

In this appendix we discuss the differences of our approach [|1^] and the treatment in [^3 
in the continuation prescription to Euclidean space with the help of a 2 x 2 dimensional 
toy modell 

The non-Hermitean structure of our Hamiltonian ( ^.5[ ) is already represented in the 
simple model: 



*This equality is not the case when e.g. matrix elements of SU(3) generators are considered |27|. 

^We thank the anonymous referee of our earlier paper fl^| for bringing this toy model to our attention. 
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Comparing the structure of h with a Dirac Hamiltonian we would like to identify the 
variable p with the "momentum" . Note that cg> 4 is real. The eigenvalues are 



ei,2 



iuj A ± sjp* - ujI (C.2) 



Then the energy functional is given by (ignoring the subtraction of the vacuum part, 
cj 4 = 0): 



eL = ^v 2 - ^ 2 r( - -, (P 2 - ^ 4 2 )/A 2 )) = (C.3) 

for the case that \p\ > |a; 4 |. Otherwise we have from eqns. ( |3.9| ) and ( |3.16| 

E* C = -^=A and = (C.4) 

i. e. a constant energy functional and thus no contribution to equation of motion. Contrary 
to this, using the prescription of ref.|24[] yields a dynamical contribution of this mode: 



* N c 



El c = j^=y/i? + <4n - ~ (P 2 + ^)/A 2 )) EL = (C.5) 



For large "momenta" p eqn. (|C.5|) obviously is the continuation of eqn ( |C.3| ) and in both 
prescriptions the contribution of states with large p is suppressed via the regularization. 
However, we see that for \p\ < |co> 4 | the two approaches are no longer identical. Especially, 
we see that in our prescription ( |C.4| ) the states with small p are not effected by the cut- 
off A in contrast to ref. |[53] . This, of course, provides a physical motivation to use our 



approach. It is also obvious from ( |C.5| ) that for \p\ < \u^\ there is still a contribution to 



the equation of motion in the formulation of ref. 0] in contradiction to the formulation 
of the determinant in Euclidean space. 
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